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We present Dalitz-plot analyses of 5+ — > K^n 7r+ and B'^ —> K^n vr" using the data sample collected at the 
T(4S') by the BABAR detector. We have found evidence for direct CP-violation in the decay — > p'^K^, with a CP- 
violation parameter A-cp = i r"^ — ~ (+44 ± 10 ± 4+?o) %, where Tg± are the decay rates. The uncertainties are 
statistical, systematic, and model-dependent, respectively. We also search for the suppressed decays B^ K~tt^tt^ 
- - . and i^+Zf +7r~ and improve upper limits on the decay branching fractions. 

oo ■ 

o ■ 
o : 

1. INTRODUCTION 
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, In the Standard Model (SM), violation of CP symmetry is a consequence of the complex phase of the Cabibbo- 
[ Kobayashi-Maskawa (CKM) quark-mixing matrix Comprehensive tests of the SM CP-violation mechanism 
require precise measurements of the three sides and three angles of the CKM unitarity triangle. Although CP- 
^ i: violation in the neutral meson [2] and meson has been well-established, it is known to be too small to 
^ account for the matter-dominated Universe, direct CP- violation would hel^^to explain the dominance of matter in 
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the Universe [J]. Direct CP-asymmetry in p^K^ has been suggested [5| and searched at BABAR and BELLE [6] 

^ through Dalitz-plot analysis of P+ K^tt^t:^ (charge conjugate decay is implied throughout this paper.) Direct 
CP-asymmetries in P ^ K*tt, KqTt, and i^T^^ are expected to be small in the SM. Moreover, the relative weak phase 
between tree and penguin diagrams in P ^ Ktttt decays is the CKM angle 7 = arg {—VudV*i,/VcdV*^). Therefore a 
J> ' set of Dalitz-plot analyses of P — > Ktttt can provide a relatively clean determination of 7 [H] . 

Compared to the penguin transitions b ~* qqd and b qqs, the decay rates for the wrong sign decays P+ ^ 
_ft'^7r+7r+ and K^K^tt^ via b dds and b ssd transitions are further suppressed by |14dl^fsl^ — O(10~^), 
[ resulting in branching fractions of ©(lO"^"*) and C(10~^^), respectively Q. Observations of these decays would be 
■ clear evidence of the b —^ dds and b — > ssd transitions. 



2. ANALYSIS TECHNIQUE 



; I ■ A number of intermediate states can contribute to 3-body P — > Ktttt decays. Their individual contributions are 
obtained from a maximum likelihood (ML) fit to the distribution of events in the Dalitz-plot formed by two invariant 
masses squared x and y of particle pairs. The total amplitudes for 3-body P and P decays are given by: 

A = A{x, y) ^ ^CjP,(x, y)- A = A{x, y) = ^CjPj(a;, y). (1) 

i j 

Where Cj is the complex coefficient for a given intermediate state j with all the weak phase dependence. The 
normalized distributions Fj describe the dynamics of the decay amplitudes and arc the product of an invariant mass 
term (relativistic Breit-Wigner in general), two Blatt-Weisskopf barrier form factors llCll. and an angular function. 
In the case of /o(980) (Kq), the mass term is replaced by the Flatte (LASS) lineshape llj. The fit fraction FF of a 
given intermediate state j with a partial decay width Tj is given by: 
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FF,^^ = -'-' ^ . (2) 
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The sum of all contributions is not necessarily unity due to the interference. The CP-asynimetry for a given inter- 
mediate state is defined as 



B^j+iSj \Cj\ +|Cj| 

The signal Dalitz-plot probability density function (PDF) is constructed as 
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Vs^gix, y,qB)= J/ : J ' , (4) 



where qs is the charge of B candidate, e = £{x,y) and e = £{x,y) are the signal reconstruction efficiencies for B 
and B events. The Dalitz-plot PDFs for continuum and BB backgrounds are constructed through two-dimensional 
histograms from simulated continuum and BB samples. 

A B candidate is reconstructed through the desired decay chains of interest. Charged kaon and pion candidates 
are identified by energy loss dE/dx information measured in a five-layer silicon vertex detector and a 40- layer drift 
chamber, and the Cherenkov angle and number of photons measured in a detector of internally reflected Cherenkov 
radiation. A tt" candidate is formed from photon pairs measured in a CsI(Tl) crystal electromagnetic calorimeter. 
The BaBar detector is described in detail elsewhere 



A B candidate is characterized by two kinematic variables: the energy-substituted mass rngs = ■\/s/4 — and 
the energy difference AE = E^ — -y/s/2, where Eg and p*^ are the center-of-mass (CM) energy and momentum of 
the B candidate, respectively, y/s is the total CM energy. Signal events peak at the nominal B mass for toes and at 
zero for AE. 

The dominant background comes from continuum production e+e^ qq, where q = u,d, s, c. This background 
is suppressed by requirements on event-shape variables calculated in the CM frame. The continuum background is 
further suppressed by exploring a neural network (NN) algorithm based on a set of kinematical variables. 

Standard extended unbinned maximum likelihood method is used to fit data. The likelihood function has the form 
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Y^rikVlix, y, TOES, A^, NN) 
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(5) 



where N = Uk is the total number of events, Uk is the fit yield of component k (k = signal, BB, and continuum). 
VI is the PDF for event i to be identified as component k. 



3. RESULTS 

3.1. Results on B+ K+n-n^ 



A Dalitz-plot analysis of the decay B^ K^Tr~TT~^ 1J| is based on a data sample of 347.5 fb~^, containing 
(383.2±4.2) X 10^ BB pairs recorded by the BaBAR detector at the T(4S') resonance at the Stanford Linear Accelerator 
Center. The Dalitz-plot variables are x = and y = to^+^_ . A phase-space non-resonant component and nine 

intermediate states K*°7:+ , K*°n+, p"K+, fo{980)K+, XcoK+, K*°n+, ujK+, /2(1270)X+, and fxK+ are included 
in the ML fit. The fit to 12,753 selected candidate events yield 4585±90±297±63 signal events and the overall direct 
CP-asymmetry of (2.8 ± 2.0 ± 2.0 ± 1.2)%, where the uncertainties are statistical, systematic, and model-dependent, 
respectively. The intermediate state ujK~^ with uj 7r"'"7r~ has noticeable effect on the p'^ lineshape and is included 
in the fit although its contribution is small. A scalar particle fx and /2(1270) are necessary to provide better fit 
to the data. The mass and width of fx are determined to be m/^ = 1479 ± 8 MeV/c^ and F/^ = 80 ± 10 MeV, 
respectively, where the uncertainties are statistical only. The fit results are summarized in Table [H The statistical 
significance of the direct CP-violation is evaluated from the differnce — 2Aln£ of the nagative log-likelihood of the 
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nominal fit and that of a fit where CP-violation parameters for the given component are set to zero, the number of 
degrees of freedom (two in this case) is taken into account. 

The total branching fraction in Table |T] is consistent with BELLE 's measurement [13|. We see evidence of direct 
CP-asymmetry of Acp ^ (+44 ± 10 ± 4+^3)% in B+ p°K+, consistent with the previous findings [6]. The 
statistical significance of the direct CP- violation effect is found to be 3.7a from the change in likelihood as described 
above. As experimental systematic uncertainties are much smaller than the statistical errors, they do not affect this 
conclusion. We have cross-checked the effect of the choice of the Dalitz-plot models on the significance. We find 
that the significance remains above Str with alternative models. The statistical significance of direct CP-violation in 
£?+ /2(1270)X+ is also above Scr, but it suffers from large model uncertainties. The direct CP-asymmetries in 
5+ ii'*'^7r+, Kq'''tt'^, and Xj^ti"^ are all consistent with the SM expectations. 



Table I: Summary of measurements of brancliing fractions (averaged over charge conjugate states) and CP asymmetries. 
Note that these results are not corrected for secondary branching fractions. The first uncertainty is statistical, the second 
is systematic, and the third represents the model dependence. The final column is the statistical significance of direct CP 
violation (DCPV) determined as described in the text. 



Mode 


Fit fraction (%) 


B{B+ Mode)(10-®) 


Acp (%) 


DCPV sig. 


iC+7r"7r+ total 




54.4 ±1.1 ±4.5 ±0.7 


2.8 ±2.0 ±2.0 ±1.2 




K*°n+; K*° K+tv- 


13.3 ± 0.7 ± 0.7 lg;g 


7.2±0.4±0.7t^-J 


+3.2±5.2± l.lti;^ 


0.9cr 


Ko°Tv+- Ko° K+Tv- 


45.0 ± 1.4 ± 1.2 lo^2® 


24.5 ± 0.9 ± 2.1 ti:? 


±3.2±3.5±2.0li:g 


1.2a 




6.54 ±0.81 ±0.58 


3.56 ± 0.45 ± 0.43 tJJ ?^ 


+44 ± 10 ±4 


3.7(7 


M980)K+; /o(980) -»7r+7r- 


18.9 ± 0.9 ± 1.7 tg;^ 


10.3 ± 0.5 ± 1.3 IJ:^ 


-10.6 ± 5.0 ± 1.1 o 


1.8a 


XcoK+; Xco tt+tt" 


1.29 ±0.19 ±0.15+0 ^3 


0.70 ± 0.10 ± 0.10 +° 02 


-14± 15±3+5 


0.5a 


K'^TT^-K^ nonresonant 


4.5 ± 0.9 ± 2.4t'i':g 


2.4±0.5± 1.3t|]j 






Kf^+. Kf K+7V- 


3.40 ± 0.75 ± 0.42 toll 


1.85 ±0.41 ±0.28 to.os 


+5±23±4l7* 


0.2a 




0.17 ±0.24 ±0.03 torn 


0.09 ± 0.13 ± 0.02 to.M 






f2(1270)K+; /2(1270) ^ n+n- 


0.91 ±0.27 ±0.11 tU'i* 


0.50 ±0.15 ±0.07^0.^9 


-85 ± 22 ± 13 If 


3.5a 


fxK^\ fx tt+tt" 


1.33 ±0.38 ±0.86 to°it 


0.73 ± 0.21 ± 0.47 tlJ JIJg 


+28 ± 26 ± 131:5 


0.6a 



3.2. Results on K+ 







A Dalitz-plot analysis of K~^tt~tt^ 15| is based on a data sample of 413 fb~^, corresponding to (454±5) x 10^ 

BB pairs produced at the T(45') resonance. The Dalitz-plot variables are x = m^+^_ and y = w^+^o . A phase-space 
non-resonant component and seven intermediate states p~{770)K'^, p~(1450)-fir+, p~{1700)K~^, K*'^Tr~ , K*^tt°, 
Kq'^tt~ , and Kq^tt^ are included in the ML fit. The fit to the data yields 4583 ± 122 signal events, where the 
uncertainty is statistical only. The results are preliminary. The decays i?*^ D ir'^ and with D K^tt^ 

and — > 7r+7r° are included in the fit as calibration modes. The fit shows that the decay i?'^ X+tt^tt" is 
dominated by KqTt and p~{770)K+. 

Isospin symmetry relates the amplitudes of /v*+7r^ and K*^^^ and B K*^tt^ and K* tt" which form 

two unitarity triangles. The orientation of the two triangles can be determined from a time-dependent Dalitz-plot 
analysis of Kgn^n^ . A Dalitz-plot analysis of i?° K^tt^tt^ can extract the magnitudes and phases of K*ti 

and K IT in B ' and B decays which are essential to determine the CKM angle 7 and will be updated soon. 



3.3. Results on B+ K-'k+'k+ and K+K+n- 

The preliminary results on the search for the suppressed decays B^ _ftr^7r+7r+ and K'^K'^n^ [l^ are based 
on a data sample of 426 fb^^ which contains (467 ± 5) x 10^ BB pairs. Since the signal yields for the two modes 
are expected to be small, no Dalitz-plot involves in the ML fit. We apply the ML fit to the 26,478 B~^ K~tt^tt~^ 
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and 7,822 S+ — > K^K^-k candidate events and find no significant signal events. We set upper limits on the decay 
branching fractions to be B{B+ K-T:+Tr+) < 9.5 x 10""^ and B{B+ -> K+K+i:-) < 1.6 x 10"'^ at 90% confidence 
level. These two upper limits have been improved by a factor of 2 and 8, respectively. 

4. SUMMARY 

We have performed a Dalitz-plot analysis of — > K~^'!t~tt~^. A scalar particle fx with rn/^ = 1479 ± 8 MeV/c^ 

and Tf-^ = 80± 10 McV and /2(1270) arc necessary to fit the data better. We find evidence for direct CP-asymmetry 
of Acp = (+44 ± 10 ± 4:tl3,)% decay B+ p"K+ with a statistical significance of 3.7a. The direct CP- 

asymmetries in iJ+ /sr*°7r+, Kq^tt^, and K^^tt^ are consistent with zero as expected. 

We have improved the Dalitz-plot analysis of K~^Tr~'K^. The determination of the magnitudes and phases 

in B^ K*Tr and B K*Tr will be updated soon which arc essential for the extraction of the CKM angle 7. 
We have improved the upper limits on the branching fractions for the SM-suppressed decays 5+ — > K~w~^n~^ and 
P+ K~^K^TT~ by a factor of 2 and 8, respectively. 
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